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ABSTRACT.  This  report  contains  the  results  of  a  parametric  study  on  the 
effectiveness  of  cluster  weapons  against  "area"  targets,  where  an  "area" 
target  consists  of  a  large  number  of  units  distributed  over  an  area 
(e.g.,  a  battalion  of  troops  in  the  field). 


The  effects  of  number  of  bomblets,  delivery  error,  ammunition  dis¬ 
persion,  and  warhead  effectiveness  are  considered.  The  results  of  the 
study  are  presented  in  handbook  form  so  that  one  may  quickly  determine 
the  effectiveness  of  any  proposed  cluster  weapon,  once  specific  values 
have  been  assigned  to  the  above  parameters. 
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INTRODUCTION 

In  an  effort  to  assist  in  the  rapid  determination  of  the  effective¬ 
ness  of  proposed  and  existing  weapons,  a  parametric  study  has  been  made 

of  such  weapon  effectiveness.  The  results  of  the  first  part  of  this 
study,  concerning  weapons  to  be  employed  against  unitary  targets,  are 
presented  in  Ref.  1.  Further  results  of  the  study,  dealing  with  the 
effectiveness  of  cluster  weapons  against  square  area  targets,  are  pre¬ 
sented  herein  in  graphical  and  tabular  form. 

The  parameters  which  were  varied  in  making  the  study  are  number  of 

bomblets  or  warheads  per  cluster  weapon,  delivery  error-*-,  ammunition  dis¬ 

persion,  radius  of  effect  of  each  warhead,  and  an  effectiveness  factor^ 
for  each  warhead. 

The  targets  under  consideration  in  this  report  are  "area"  targets 
as  opposed  to  "unitary"  targets.  A  "unitary"  target  is  a  target  such  as 
a  tank  or  a  ship  which  may  be  regarded  as  a  unit  for  purposes  of  tactical 
warfare.  The  "area"  target  consists  of  a  large  number  of  units  distrib¬ 
uted  over  an  area.  The  results  presented  in  this  report  are  based  on 
the  assumption  that  this  area  is  square  in  form. 

The  results  of  studies  on  problems  similar  to  the  one  treated  here 
are  found  in  Ref.  2  and  3.  The  basic  assumptions  found  in  Ref.  2  are 
significantly  different  from  the  basic  assumptions  of  the  study  reported 
herein;  hence,  the  results  of  the  present  study  and  those  of  Ref.  2  are 
not  directly  comparable.  In  Ref.  3j  the  basic  assumptions  are  such  that 
a  comparison  may  be  made  with  some  of  the  results  of  the  present  study. 
The  comparison  appears  in  the  third  section  of  this  report. 


■METHOD  OF  CALCULATION 

It  is  desired  to  determine  the  fraction  T  of  a  two-dimensional  tar¬ 
get  which  may  be  expected  to  be  destroyed  by  a  cluster  of  bomblets.  The 
following  assumptions  are  made: 


1  The  expression  "delivery  error"  as  used  in  this  report  includes 
all  errors  in  the  system  resulting  in  the  displacement  of  the  projected 
impact  point  of  the  unseparated  cluster  from  its  aimpoint.  Specifically, 
the  effects  of  aim  or  fire-control  error  and  the  dispersion  of  the  clus¬ 
ter  weapon  as  a  whole  are  included  In  the  delivery  error. 

2  The  meaning  of  the  tenra  "effectiveness  factor"  as  used  in  this  re¬ 
port  is  explained  in  the  section  on  the  method  of  calculation. 
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1.  The  target  is  a  square  with  its  sides  parallel  to  the  axes  of  a 
Cartesian  coordinate  system. 

2.  The  projected  impact  point  of  the  cluster  is  normally  distributed 
about  its  aimpoint,  due  to  delivery  error. 

3.  The  bomblets  are  normally  distributed  about  the  projected  im¬ 
pact  point,  due  to  ammunition  dispersion. 

4.  The  delivery  error  is  independent  of  the  ammunition  dispersion, 

5.  Each  borablet  is  uniformly  effective  throughout  a  square  region 
centered  on  the  impact  point  of  the  borablet. 

Since  it  is  desired  to  take  account  of  cumulative  damage  which  might 
result  from  several  bomblets  impacting  in  proximity  to  one  another,  a 
Monte  Carlo  method  of  evaluation  is  used.  The  method  involves  the  use 
of  random  normal  numbers  in  repeated  evaluations  of  the  weapon,  the  pro¬ 
cess  being  continued  until  the  expected  error  in  the  estimate  of  weapon 
effectiveness  is  satisfactorily  small. 

The  target  is  defined  so  that  it  lies  within  a  "playing  field"  which 
is  a  100- X  100-unit  square.  The  location  of  any  lattice  point  in  this 
square  is  specified  by  an  ordered  pair  of  positive  integers  which  repre¬ 
sent  the  coordinates  of  the  point  in  a  Cartesian  coordinate  system. 

Four  quantities  are  required  to  characterize  the  target:  the  co¬ 
ordinates  of  the  lower  left  corner,  the  length  of  the  sides,  and  a  unit 
value.  For  instance,  if  the  sides  of  the  target  are  3  units  in  length 
with  a  unit  value  of  10,  each  square  unit  in  the  target  will  have  a 
"value"  of  10  so  that  the  total  "value"  of  the  target  is  90*  It  is  con¬ 
sidered  that  each  square  unit  of  a  target  is  centered  on  a  lattice  point 
of  the  "playing  field"  and  its  value  is  concentrated  there.  That  is,  if 
a  target  is  defined  with  lower  left  corner  at  (2,  6)  and  lengths  of  sides 
one  unit  each,  only  the  point  (2,  6)  of  the  "playing  field"  will  be  con¬ 
sidered  to  be  occupied  by  the  target. 

If  a  lattice  point  of  the  "playing  field"  is  included  in  the  target, 
the  point  is  said  to  be  a  target  point  and  the  unit  value  of  the  target 
is  assigned  to  the  point.  The  total  value  of  the  target  is  seen  to  be 
the  sum  of  the  values  of  the  target  points.  The  fraction  T  of  this  total 
value  which  is  expected  to  be  destroyed  by  the  cluster  weapon  is  determined 
by  the  calculation  procedure. 

Two  quantities  are  required  to  characterize  the  effectiveness  of  each 
bomblet:  the  effectiveness  distance  R  and  the  effectiveness  factor  E. 

The  values  of  these  quantities  are  the  same  for  all  bomblets  in  a  given 
weapon.  The  square  area  of  effectiveness  of  a  bomblet  (which  will  hence¬ 
forth  be  named  the  impact  area  of  the  bomblet)  includes  those  lattice 
points  which  lie  within  an  x-distance  of  R  from  the  impact  point  and 
which  also  lie  within  a  y-distance  of  R  from  the  impact  point.  The  im¬ 
pact  point  will  always  be  a  lattice  point.  For  instance,  if  R=l,  the 
impact  area  of  each  borablet  consists  of  the  nine  lattice  points  centered 
on  the  impact  point  of  that  bomblet  (three  rows  of  three  points  each). 
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The  effectiveness  of  an  individual  borablet  is  accounted  for  in  the 
following  way.  If  a  point  of  the  impact  area  of  a  borablet  coincides 
with  a  target  point,  the  value  of  the  target  point  is  multiplied  by  the 
effectiveness  factor  E  (E  £  l),  the  product  is  added  to  a  tally  of  tar¬ 
get  value  destroyed,  and  the  value  of  the  target  point  is  reduced  by  the 
value  of  the  product.  If  the  resulting  value  is  less  than  0.5,  the  tar¬ 
get  point  is  given  the  value  zero  (at  which  time  it  loses  its  status  as 
a  target  point). 

Let  us  assume  that  a  target  has  been  defined  and  the  characteristics 
of  the  borablets  (impact  area,  effectiveness,  dispersion),  as  well  as  the 
number  of  bomblets  in  the  weapon,  the  aim  point  of  the  weapon^,  and  the 
delivery  errors  have  been  determined.  The  effectiveness  of  the  weapon 
against  the  target  is  determined  as  follows. 


Two  random  normal  numbers,  r  and  r  ,  are  selected  from  a  list  of 
random  normal  numbers  with  unit  standard^deviation  and  zero  mean4.  Then 
the  coordinates  ax  and  ay  of  the  actual  aimpoint  for  the  weapon  are  given 
by 


ax  *  Ax  +  rl*  F. 


X 


ay  -  Ay  +  Tp  •  <T~1 :>  , 


F, 


(1) 


where  Ax  and  Ay  are  the  coordinates  of  the  intended  aimpoint  of  the  wea¬ 
pon,  and  CT px  and  C p  are  the  standard  deviations  of  the  delivery 
error  in  the  x-  and  y-dlrections.  Note  that  ax  and  a y  are  not  necessarily 
integers. 


Now  two  more  random  normal  numbers,  r^  and  rt ,  are  generated,  and 
the  coordinates  Iq^,  and  Iq  of  the  actual  impact  point  of  the  first  bomb- 
let  are  given  by:  y 

^lx  =  (ax  +  r3*0"R)  rounded 

I ly  =  (ay  +  r4-cnty.}  rounded, 

where  Ilx  and  Iq  have  been  rounded  to  integers;  ax  and  ay  are  the 

quantities  determined  by  Eq.  1,  and  CT  px  and  cr  Ry  are  the  standard 
deviations  of  the  ammunition  (bomblets)  in  the  x-  and  y-directions. 


3  For  the  calculations  of  this  report,  the  aimpoint  of  the  weapon 
was  chosen  to  be  the  center  of  the  target. 

If  this  method  is  pi-ograramed  for  a  digital  computer,  subroutines 
are  available  which  generate  random  normal  numbers  with  the  desired 
characteristics . 
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Once  Ii  and  It.  have  been  determined,  the  impact  area  for  the  first 
bomblet  is  compared  with  the  target  area  to  determine  whether  the  impact 
area  of  this  bomblet  includes  any  target  points.  If  any  target  points 
are  so  included,  the  value  destroyed  is  accounted  for  as  described  pre¬ 
viously  and  the  target  is  modified  by  appropriately  reducing  the  value 
of  these  target  points. 

The  actual  impact  point  of  the  second  bomblet  is  now  determined  in 
a  similar  manner,  account  is  taken  of  any  damage  the  second  bomblet  does 
to  the  target,  and  the  process  is  repeated  until  all  bomblets  have  been 
accounted  for.  For  the  ith  bomblet,  equations  analogous  to  Eq.  2  yield 

Iix  -  (ax  +  r  j  *  mounded 

*ly  *  (ay  +  rk.  fffcy)  rounded, 

where  r  j  and  rk  are  two  random  normal  numbers.  Note  that  for  each  bomb- 
let  theJsame  values  of  ax  and  ay  are  used. 

When  a.11  bomblets  have  been  accounted  for,  the  accumulated  value 
destroyed  is  divided  by  the  original  total  target  value  and  the  result 
is  the  first  element  in  the  set  of  evaluations  of  the  weapon  against  the 
target.  The  target  is  now  redefined  in  its  original  state  and  the  pro¬ 
cess  begins  again  with  the  selection  of  another  actual  aimpoint  for  the 
weapon  (Eq.  l).  After  many  repetitious  of  the  procedure,  the  final  an¬ 
swer  T  is  taken  as  the  average  of  the  results  of  the  several  trials. 


In  order  to  determine  the  point  at  which  to  stop  the  process  and 
give  an  estimate  of  the  error  in  the  final  result,  the  quantity 
where  a  is  the  average  value  of  the  elements  of  {  a  }  ,  is  determined  for 
set  {  a  }  of  results  of  the  trials.  When  (J\  becomes  less  than  a  speci¬ 
fied  value,  the  process  is  stopped  and  we  then  have  a  probability  of  a- 
bout  0.95  that  the  correct  result  T  satisfies  |T  -  a|  <  2 It  is 
known  that 


cr-  = 


Ok 


yiET 


V-Zifai-S)1 


n 


(M 


Equation  (4)  may  be  used  to  evaluate  C7g. 


If  the  number  of  bomblets  per  cluster  weapon  is  at  all  sizable,  it 
is  apparent  that  the  above  procedure  is  too  tedious  for  hand  computation. 
However,  the  problem  is  suitable  for  programming  on  a  high-speed  digital 
computer,  and  this  has  been  done:  an  IBM  709  program  to  carry  out  the 
computation  procedure  has  been  written.  A  complete  description  of  this 
program  will  be  found  in  Ref.  4.  The  program  was  used  to  compute  the 
estimates  of  weapon  effectiveness  which  are  presented  in  this  report. 
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COMPARISON  WITH  ANOTHER  METHOD 

As  mentioned  earlier,  some  work  of  a  similar  nature  to  that  presented 
in  this  report  is  described  in  Ref.  3*  The  assumptions  of  Ref.  3  are 
somewhat  different  from  those  of  this  report,  and  the  method  used  there 
is  analytic  in  nature  rather  than  being  a  Monte  Carlo  method.  However, 
it  is  possible  to  compare  some  results  of  the  two  methods,  and  this 
section  contains  the  comparison. 

The  principal  difference  between  the  assumptions  of  this  report  and 
those  of  Ref.  3  is  that  no  differentiation  is  made  in  Ref.  3  between 
delivery  errors  and  ammunition  dispersion.  The  average  fraction  of  a 
square  target  which  is  covered  at  least  once  by  up  to  10  individually- 
dropped  bombs  (each  having  a  square  impact  area)  Is  calculated  in  Ref.  3* 
Since  these  bombs  are  not  dropped  in  a  cluster,  there  is  no  need  in 
Ref,  3  to  assume  a  separate  delivery  error. 

The  results  of  this  report  which  are  comparable  to  data  in  Ref.  3 
arc  those  for  which  CHp  =  0,  N  *  10  and  E  »  1.  It  is  assumed  both  in  this 
report  and  in  Ref.  3  that  O^Rx  =  CT" Table  1  contains  estimates  of 
weapon  effectiveness  from  both  Ref.  S^and  this  report.  The  impact  areas 
of  the  bomblets  range  from  the  target  area  to  .0156  (one  sixty-fourth) 
of  the  target  area.  The  quantity  AE  in  Table  1  is  the  impact  area  of 
the  borablet  divided  by  the  target  area.  The  quantity  O'"  r  is  expressed 
in  units  of  target  widths.  It  is  the  common  value  of  CTrx  and  CTpy. 

The  somewhat  unusual  values  of  O'r  arise  as  a  result  of  the  method  of 
tabulation  of  the  results  in  Ref.  3. 

TABLE  1.  Expected  Fraction  of  Target  Destroyed  10  Bomblets 


r  .0156 

,  .0625 

.25 

1.0 

0 

.0156 

.0625 

.25 

1.0 

.053 

,0628a 

.062 

.1509 

.163 

.4201 

.445 

.9999 
•  1.0 

.106 

.1059 

.102 

.2511 

.256 

.5972 

.64o 

.9997 

1.0 

.212 

.1329 

.131 

.3930 

.377 

.7960 

.807 

,  .9990 

1.0 

.425 

.0863 

.091 

.3007 

.309 

.7494 

.756 

.9937 

.998 

.850 

.0303 

.1158 

.3884 

.410 

.8583 

.879 

a  The  upper  number  of  each  pair  is  from  Ref.  3  and  the  lower  number 
is  from  this  report. 
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It  is  seen  from  an  inspection  of  Table  1  that  the  agreement  between 
the  two  methods  Is  quite  good.  In  most  cases  the  difference  between  the 
results  is  less  than  0.01,  and  the  few  cases  where  the  difference  is 
greater  are  partly  due  to  artistic  license  in  drawing  the  curves  of  this 
report.  The  maximum  discrepancy  is  0.04. 


ASSUMPTIONS  AND  PROGRAM  DESCRIPTION 

In  making  the  calculations  for  this  report,  several  assumptions  were 
made  in  addition  to  those  listed  in  the  description  of  the  method  of  cal¬ 
culation.  For  completeness,  all  major  assumptions  and  conditions  are 
listed  below.  These  assumptions  were  made: 

1.  The  target  is  a  square  with  its  sides  parallel  to  the  axes  of  a 
Cartesian  coordinate  system. 

2.  The  aimpoint  of  the  weapon  is  the  center  of  the  target. 

3.  The  projected  impact  point  of  the  cluster  is  normally  distributed 
about  its  aimpoint,  due  to  delivery  error. 

4.  The  bomblets  are  normally  distributed  about  the  projected  impact 
point,  due  to  ammunition  dispersion. 

5.  The  delivery  error  is  independent  of  the  ammunition  dispersion. 

6.  Each  bomblet  is  uniformly  effective  throughout  a  square  region 
centered  on  the  impact  point  of  the  bomblet. 

7.  ffFv  -  °Ry  _  i 

Qfx  “  ^ 

8.  All  target  points  have  an  initial  unit  value  of  10. 

Since  assumption  7  indicates  that  Of  y  =  Of  x  and  CTRy  =  <?“Rx»  it 
will  not  be  necessary  to  distinguish  between  the  two  directions  when 
delivery  errors  and  ammunition  dispersion  are  discussed.  Therefore,  the 
symbols  CTp  and  CTr  will  henceforth  be  used  when  reference  is  made  to 
delivery  errors  and  ammunition  dispersion. 

Assumption  8  increases  the  sensitivity  of  the  method  for  cases  in 
which  E<  1. 

The  following  parameters  were  varied  in  making  this  study:  Of  ,  C7”b, 
N,  R,  and  E.  The  unit  for  <7"p,  and  R  is  target  length. 

Due  to  the  way  in  which  R  is  defined,  the  area  of  effectiveness  of 
a  bomblet,  Ag,  is  given  by 

A£  =  (2R  +  if.  (5) 

In  making  the  calculations  for  this  report,  values  of  R  were  chosen  such 
that  Ag  assumed  values  of  one  sixty- fourth,  one  sixteenth,  one  fourth, 
and  one  target  area.  The  decimal  expansions  of  these  values  appear  in 
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Table  2  below. 

Table  2  contains  the  values  of  &F,  N,  Ag,  and  E  (bomblet  effect¬ 
iveness  factor)  which  were  used  in  the  calculations. 

TABLE  2.  Values  of  <Jp  ,  N,  AEj  and  E  Used  in  the  Study 

Parameter  Values 

E  0.5,  1 

N  1,  10,  50,  100 

C7p  (Target  lengths)  0,  0.25,  0.5,  1,  2,  4 
Ag  (Target  areas)  0.015,  0.06,  0.25,  1 

For  calculations  involving  a  given  value ^  of  Of,  several  values  of 
cTg  were  chosen  in  the  range  0<-  Of  <  3  F*  In  this  way  the  variation 
of  T  with  was  determined.  It  was  assumed  that  the  ammunition  dis¬ 
persion  is  easier  to  control  than  the  delivery  error,  and  when  optimized 
values  of  T  are  discussed,  the  optimization  is  with  respect  to  <T r.  That 
is ,  when  the  data  are  crossplotted  in  Appendixes  C  and  D,  the  maximum 
value  of  T  from  each  curve  in  Appendix  A  is  used  for  the  crossplotting. 


USE  OF  THE  DATA 


Appendix  A  contains  the  basic  data  of  the  study.  The  results  for 
N=1  were  obtained  by  an  analytical  method  and  are  tabulated.  The  re¬ 
sults  for  the  other  values  of  N  are  presented  in  graphical  form.  Each 
graph  of  Appendix  A  contains  several  curves  of  fraction  of  target  de¬ 
stroyed  T  vs.  ammunition  dispersion  Of.  The  various  curves  were  ob¬ 
tained  by  use  of  different  values  of  Ag  and  E  for  fixed  values  of  CTp  and 
N.  Results  are  not  given  for  all  values  of  Ag,  since  in  some  cases  the 
delivery  error  was  so  large  and/or  N  was  so  small  that  the  smaller  values 
of  Ae  gave  negligible  results. 

Using  the  method  of  optimization  outlined  in  the  last  section  (i.e., 
choosing  a  value  of  such  that  T  is  maximum  for  a  given  curve),  the 

maximum  values  of  fraction  of  target  destroyed  are  tabulated  in  Ap¬ 
pendix  B.  These  maximum  values  of  T  are  plotted  against  N  in  Appendix  C 
to  show  the  effect  of  number  of  bomblets,  and  they  are  plotted  against 
CTf  In  Appendix  D  to  show  the  effect  of  delivery  error. 


Elementary  Examples 


Since  several  parameters  are  involved  in  the  calculations,  it  is 
clear  that  the  simplest  evaluations  to  make  will  be  those  for  which  the 
parameters  of  interest  coincide  with  values  of  those  parameters  used  in 


5  For  the  calculations  in  which  (7"  =0,  the  values  of  Cff  were 

chosen  in  the  range  0<  CT^  <  2  target  lengths. 


7 


NAVWEPS  REPORT  7641 


the  calculations*  For  instance,  if  one  is  interested  in  determining  T 
for  a  weapon  containing  50  bomblets  against  a  100-x  100-ft  target  where 
<7"p»50  ft,  C~ r=50  ft,  Ag=,25,  and  E=.5j  he  proceeds  as  follows.  Since 
the  target  is  100  ft  wide  and  the  units  for  p  and-  <T~ p  are  target 
widths  W,  <T~ p=*5W,  and  CT rs*5W  in  this  case.  The  data  for  N=50  and 
(JT p=.5W  are  found  in  Fig.  A9.  Using  the  dashed  curves  of  this  figure 
for  E».5  and  the  second  from  the  highest  of  these  (keyed  with  a  square) 
for  Ae=.25>  we  find  that  T=G,70  when  <J~ r=.5W.  Incidentally,  it  happens 
that  O'] r=.5W  gives  the  greatest  value  of  T  for  this  particular  curve. 

In  other  words,  5W  is  the  optimum  value  of  C~ 3  for  N»50#  C7f=.5W, 

Ap=.25,and  E=.5. 

If  it  can  be  assumed  that  the  optimum  value  of  0~\ r  can  be  achieved 
for  a  given  weapon,  it  Is  no  more  difficult  to  evaluate  T  when  either  N 
or  CTp  (but  not  both)  deviates  from  a  value  for  which  calculations  were 
made.  In  this  case.  Appendix  C  or  D  may  be  used.  For  example,  if  N=4q 
in  the  previous  example  and  the  values  of  the  other  parameters  are  the 
same,  we  use  Fig.  C3  (  CT p=.5W)  to  get  T0=0,65  from  the  dashed  curve  keyed 
with  a  square6.  If  ft  (.75*0  and  N=*50,  Fig.  D3  (N=50)  gives 

To=0.55. 


Interpolation 

In  the  previous  examples,  not  much  more  complexity  is  introduced  by 
allowing  E  to  assume  an  intermediate  value  and  using  the  following  para¬ 
bolic  interpolation  formula: 

(4m  -  n)E  +  (2r,  -  4m)E2,  (6) 

where  m  and  n  represent  values  of  T  when  E=0.5  and  1,  respectively,  and 
the  values  of  all  other  parameters  are  fixed.  It  is  seen  from  Eq.  6 
that  T=m  when  E=0. 5  and  T=n  when  E=l,  so  that  the  curve  whose  graph  is 
Eq.  6  passes  through  the  points  for  which  calculations  have  been  made?. 

To  illustrate  the  use  of  Eq.  6,  let  us  consider  the  second  of  the 
above  examples.  We  found  that  To=0.65  when  N=4o,  <7p=0.5W,  Ae=.25, and 
E=.5.  Let  us  find  T0  when  E=,7  and  the  values  of  the  other  parameters 
are  the  same  (except  that  the  optimum  value  of  CT"r  may  be  slightly  dif¬ 
ferent  due  to  the  change  in  E).  .  We  consult  Fig.  C3  again  and  find  that 
To=0.82  when  E=1  (solid  curve  keyed  with  a  square).  Hence,  we  have 
m=0.65  and  n=0.82  for  this  application  of  Eq.  6,  and  our  interpolation 


6  The  notation  T0  denotes  optimization  with  respect  to  cTr. 

?  It  has,  of  course,  not  been  established  that  Eq.  6  represents 
the  proper  curve  for  interpolation.  However,  use  of  this  curve  appears 
to  the  author  to  yield  more  realistic  results  than  would  be  the  case  if 
linear  interpolation  were  used. 
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formula  for  this  example  is 

Tq=  [  4(0.65)  -  O.82]  E  +  [2(0.82)  -  4 (0.65 )_7  E2 
3  I.78E  -  O.96E2.  (7) 

Substituting  E=0.7  into  Eq.  7>  we  find  that  To=0,78  when  E=0.7» 

If  it  is  desired  to  interpolate  on  area  of  effectiveness  Ae.,  the 
situation  becomes  more  complicated.  Perhaps  the  best  procedure  would  be 
to  evaluate  T  (or  T0)  for  all  values  of  Ag  available  in  this  report  and 
sketch  a  graph  of  T  (or  T0)  versus  Ag.  Then  the  interpolation  can  be 
made  by  reading  the  graph  for  the  appropriate  value  of  Ag. 

Graphical  interpolation  is  somewhat  laborious.  If  a  numerical 
method  of  interpolation  on  Ag  is  desired,  perhaps  the  simplest  and  best 
is  linear  interpolation.  The  principal  region  in  which  such  a  method  is 
Inaccurate  is  the  region  between  Ag=0.25  and  Ag=l,  and  any  errors  in  this 
region  will  be  on  the  conservative  side.  Therefore,  if  the  results  of 
this  report  are  to  be  used  for  preliminary  estimates  of  weapon  effective¬ 
ness,  linear  interpolation  on  Ag  should  not  introduce  catastrophic  errors. 

Let  us  compare  the  results  of  a  graphical  interpolation  on  Ag  with 
the  results  of  a  linear  interpolation  on  the  same  parameter.  We  shall 
take  E»l,  N=4o,  <T jr«0.5W,  <T r  its  optimum  value,  and  we  shall  find  Tc 
for  AgaO.5.  Prom  Fig.  C3  we  find  that  To=0.l8  when  Ae*0.015>  To=0.46 
when  Ajr=0*062,  To=0.82  when  Ag=0.25,  and  To=0.99  when  Ae=1.  If  we  plot 
these  points,  we  find  that  Tq=0.94  when  AgaO.5.  If  we  use  linear  inter¬ 
polation  on  Ag  we  get  TqjsO.88  when  Ae=0.5.  If  instead  of  doing  a  linear 
interpolation  on  Ag  we  do  it  on  '/Ag  (the  length  of  a  side  of  the  square 
area  of  effectiveness),  we  get  To=0.89  when  -yAg=  -/o. 5  =  0.7*  This 
latter  result  is  in  slightly  better  agreement  with  the  result  obtained 
graphically. 

If  intermediate  values  must  be  used  for  all  parameters,  the  inter¬ 
polation  problem  becomes  quite  formidable.  We  conclude  this  section  of 
the  report  with  example  in  which  we  evaluate  T0  for  E=0.7j  Ag=0*5>  N=40, 
and  CTj«0.75.  One  way  of  proceeding  is  as  follows.  Taking  E=0.5*  we 
evaluate  T0  for  the  four  values  of  Ag  of  this  report  when  N=40  and  (7p=* 
0.5W.  We  use  Fig.  C3  and  get: 


Ak 

E 

Gj? 

N 

0.11 

0.015 

0.5 

0.5W 

4o 

0.28 

0.062 

0.5 

0.5W 

40 

0.65 

0.25 

0.5 

0.5W 

4o 

0.90 

1.00 

0.5 

0.5W 

40 

We  now  repeat  the  process  for  CJp=W  and  the  same  values  for  the 
other  parameters  (Fig.  C4)j 
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Ae 

E 

Op 

N 

0.04 

0.015 

0.5 

w 

4o 

0.13 

0.062 

0.5 

w 

40 

0.37 

0.25 

O'.  5 

w 

40 

0.69 

1.00 

0.5 

w 

40 

If  we  examine  the  curves  of  Fig.  D3  (for  N=50),  we  see  that  linear 
approximations  to  the  curves  are  satisfactory  in  the  range  0.5  ^  l. 

Therefore  we  perform  linear  interpolations  on  the  data  tabulated  above 
to  evaluate  T0  at  Of  -0.75W: 


Za 

ae 

E 

Of 

N 

0.08 

0.015 

0.5 

0.75W 

40 

0.21 

0.062 

0.5 

0.75W 

40 

0.51 

0.25 

0.5 

0.75W 

4o 

0.80 

1.00 

0.5 

0.75W 

40 

Wg  are  now  ready  to  plot  T0  vs.  Ag  to  evaluate  TQ  when  AgaO.5.  The 
result®  is  T0  =*0.67. 


In  order  to  interpolate  on  E,  it  will  be  necessary  to  begin  again 
U3ing  E=sl. 


By  use  of 

Fig.  C3 

we  have 

£0 

^E 

E 

Of 

w 

0.18 

0.015 

1.0 

0.5W 

40 

0.46 

0.062 

1.0 

0.5W 

4o 

0.82 

0.25 

1.0 

0.5W 

40 

0.99 

1.00 

1.0 

0.5W 

4o 

get  from 

Fig.  C4 

To 

Ag 

E 

Of 

N 

0.09 

0.015 

1.0 

w 

40 

0.22 

0.062 

1.0 

w 

40 

0.48 

0.25 

1.0 

w 

4o 

0.87 

1.00 

1.0 

w 

4o 

We  again  do  linear  interpolations  on  CT p  to  obtain 


T 

ia 

ae 

E 

Op 

N 

0.14 

0.015 

1.0 

0.75W 

4o 

0.34 

0.062 

1.0 

0.75W 

4o 

0.65 

0.25 

1.0 

0.75W 

4o 

0.93 

1.00 

1.0 

0.75W 

40 

®  If  linear  interpolation  on  ifkp  were  used  instead  of  plotting  the 

eraph,  it  would  not  have  been  necessary  to  evaluate  T0  for  Ae=0.01  and  Ae*= 
.06.  The  result  of  the  linear  interpolation  would  be  To=0,63  when 

-y^E=o.7. 
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Another  plot  is  made  of  T0  vs  Ae  to  evaluate  T0  when  Ae=>0.  5.  This  time 
the  result  is  Toss0.8l, 

We  are  now  ready  to  interpolate  on  E,  using  Eq.  6.  Having  To**0,67 
when  E**0. 5  and  To=0.8l  when  E=«1.0j  we  see  that  m=0*67  and  n**0.8l.  There¬ 
fore  Eq.  6  yields 

T„=n+(0.67)  -  0.8lJ  E  +  [2(0.81)  -  4  (0.67)1  E2 

2  (8) 

»  I.87E  -  1.06E  . 

Substituting  E=0.7  into  Eq.  8,  we  reach  the  final  result  of  TQ«0.78 
when  E=0.7>  AjpO.5,  N=40,  and  CTp  =0.75*  By  a  separate  calculation,  the 
correct  answer  for  this  example  was  found  to  be  To=0.75* 


CONCLUSION 

The  results  of  a  parametric  study  to  determine  the  effectiveness  of 
cluster  weapons  against  square  area  targets  have  been  presented  and  dis¬ 
cussed  in  this  report.  The  method  of  calculation  of  the  fraction  of  tar* 
get  destroyed  has  been  described  and  examples  have  been  given  to  illus¬ 
trate  the  use  of  the  material.  The  applicability  of  the  data  presented 
in  this  report  is  quite  general,  since  the  quantities  involved  in  the 
computations  are  dimensionless  and  independent  of  any  assumptions  re¬ 
garding  specific  targets  or  warheads. 
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Appendix  A 
BASIC  DATA 

TABLE  Al.  Expected  Fraction  of  Target  Covered  by  One  Bomb 


1.0 

- 1 

0,25 

0.0625 

0.0156 

0.25 

o.64oo 

0.2100 

0.0569 

0.0156 

0.50 

0.3700 

0.1075 

0.0288 

0.0072 

1.00 

0.1350 

0.0360 

0.0090 

0.0023 

2.00 

0.0380 

0.0098 

0.0024 

0.0006 

4.00 

0.0100 

0.0025 

0.0006 

0.0002 

It  should  be  noted  that  the  quantity  tab-  lated  above  is  fraction  of 
target  covered.  If  it  is  desired  to  evaluate  the  fraction  of  target  de¬ 
stroyed  by  a  particular  bomb,  the  numbers  in  the  table  should  be  multi¬ 
plied  by  the  effectiveness  factor  E  of  the  bomb. 

Since  only  one  bomb  is  involved,  the  delivery  error  <T p  and  the 
ammunition  dispersion  CTr  need  not  be  distinguished,  and  an  overall 
error  CT  will  suffice.  The  relation 

C72  -  of*  o£ 

is  the  proper  way  to  combine  CT p  and  CT r  in  this  case. 
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PIG.  Al.  fraction  of  Target  Destroyed  vs.  Ammunition  Dispersion 
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FIG.  A5.  Fraction  of  Target  Destroyed  vs.  Ammunition  Dispersion. 
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FIG.  A6.  Fraction  of  Target  Destroyed  vs.  Ammunition  Dispersion. 
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Fraction  of  Target  Destroyed  vs.  Ammunition  Dispersion 
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FIG.  A8e  Fraction  of  Target  Destroyed  vs.  Ammunition  Dispersion, 

CTV  =  0.25W,  N  “50 
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FIG,  A15.  Fraction  of  Target  Destroyed  vs.  Ammunition  Dispersion 

cr_  =  0.5W,  N-100 
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FIG.  A18.  Fraction  of  Target  Destroyed  vs.  Ammunition  Dispersion. 
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Appendix  B 

OPTIMISATION  OF  BOMBLET  DISPERSION 
TABLE  Bl.  Optimum  Values  of  Fraction  of  Target  Destroyed 

EaO.5 


Ae  1.00 

N=10 

0.25 

0.0625 

0,0156 

0 

1.00 

0.62 

0.24 

0.07 

0.25 

0.86* 

0.48 

0.19 

0.06 

0.50 

0.68 

0.32 

0.10 

0,03 

1.00 

0.35 

0,l4 

0.03 

0.01 

2.00 

0,l4 

0.04 

0.02* 

0,00* 

A. 00 

o.o4 

0.01* 

0.00* 

0.00* 

N»50 

0 

1.00 

0.95 

0.66 

0.28 

0.25 

0.97* 

0.87 

0.52 

0,22 

0.50 

0.93 

0.70 

0.33 

0,l4 

1.00 

0.?4 

0.42 

0.16 

0.05 

2.00 

0.43 

0.15 

0.05 

0.02 

4.00 

0.17 

0.05 

0.01* 

0,00* 

NaslOO 


0 

1.00 

1.00 

0.82 

0.46 

0.25 

0.99* 

0.95 

0.71 

0.35 

0.50 

0.98  . 

0.84 

0.53 

0,21 

1.00 

0.88 

0.57 

0.26 

0,08 

2.00 

0.60 

0.25 

0.09 

0.02 

4.00 

0.27 

0.10 

0.03 

0.01* 

Values  marked  by  an  asterisk  00  were  determined  by  interpolation 
or  extrapolation.  The  quantity  UF  is  expressed  in  units  of  target 
widths.  The  "optimum  values  of  fraction  of  target  destroyed"  are  the 
values  of  T  for  which  is  optimized  as  described  in  the  report. 
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TABLE  B2.  Optimum  Values  of  Fraction  of  Target  Destroyed 


E=1.0 

N=10 


AE  1.00 

0.25 

0.0625 

0.0156 

0 

1,00 

0.83 

0.40 

0.13 

0.25 

0.93* 

0.69 

0.31 

0.11 

0.50 

0.83 

0.47 

0.17 

0.06 

1.00 

0.52 

0.21 

0.07 

0.03 

2.00 

0.23 

0.10 

0.02 

0.01 

4.00 

0.09 

0.03 

0.01* 

0.00* 

N=50 


0 

1.00 

1.00 

0.84 

0.46 

0.25 

1.00* 

0.98 

0.72 

0.36 

0.50 

1.00 

0.86 

0.51 

0.22 

1.00 

0.91 

0.55 

0.25 

0.10 

2.00 

0.60 

0.27 

0.08 

0.03 

4.00 

0.29 

0.09 

0.03 

0.01* 

N»100 

0 

1.00 

1.00 

0.96 

0.65 

0.25 

1.00* 

1.00 

0,87 

0.52 

0.50 

1.00 

0.95 

0.68 

0.33 

1.00 

0.97 

0.71 

0.39 

0.15 

2.00 

0.76 

0.4l 

0.16 

0.06 

4.00 

0.40 

0.16 

0.05 

0.02* 

Values  marked  by  an  asterisk  (*)  were  determined  by  interpolation 
or  extrapolation.  The  quantity  (tf  is  expressed  in  units  of  target 
widths . 
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Appendix  C 

EFFECT  OF  NUMBER  OF  BOMBLETS 
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PIG.  C2.  Optimum  Fraction  Destroyed  vs.  Number  of  Bomblets. 

cri=  0.25W 
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FIG.  C4.  Optimum  Fraction  Destroyed  vs.  Number  of  Bamblets. 
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FIG.  G5.  Optimum  faction  Destroyed  vs.  Number  of  Bomblets 
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FIG.  C6.  Optimum  Fraction  Destroyed  vs.  Number  of  Bomblets 
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Appendix  D 

EFFECT  OF  DELIVERY  ERROR 
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PIG*  Dl*  Optimum  Fraction  Destroyed  vs*  Delivery  Error* 
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FIG.  D2*  Optimum  Fraction  Destroyed  vs.  Delivery  Error 
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FIG.  D3-  Optimum  Fraction  Destroyed  vs.  Delivery  Error 
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FIG,  D4.  Optimum  Fraction  Destroyed  vs.  Delivery  Error 
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